Charge transfer statistics and entanglement in normal-quantum dot-superconductor 

hybrid structures 
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We analyze the full counting statistics (FCS) of a single-site quantum dot coupled to multiple 
metallic electrodes in the normal state and a superconductor for arbitrary transmission. We present 
an analytical solution of the problem taking into account the full energy dependence of the transmis- 
sion coefficient. We identify two transport processes as sources of entanglement between the current 
carriers by observing positive cross current correlations. Furthermore, we consider ferromagnetic 
electrodes and show how they can be used as detectors in experiments violating the Bell-Clauser- 
Horne-Shimony-Holt inequality. 
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Correlations represent one of the main ingredients 
of quantum mechanics. They may persist even if the 
particles are causally disconnected which allows for ex- 
perimental Bell inequality tests, quantum cryptography 
and quantum teleportationi. In quantum optics using 
e. g. parametric downconversion 2 many interesting 
experiments investigating the fundamentals of quantum 
mechanics became possible^. Solid state entanglers 
have only recently been shown to work experimentally 
using a superconducting finger coupled to normal leads 
directly^ - - or via two quantum dots^. They rely on 
the non-local or crossed Andreev reflection (CAR) in 
multi-terminal structures with a superconducting lead 
where an incident electron from one of the electrodes is 
reflected by the superconductor into another electrode. 
This process produces entanglement between the current 
carriers in the two leads since the correlations of the spin 
singlet (Cooper pair) are transferred to the spatially 
separated electrons. The nature of CAR has been 
investigated further using ferromagnetic contacts 9 -*^. 
Recently, the nonlocal conductance in superconductor 
hybrid structures has been extensively studied theoreti- 
cally using different approaches^— . On the one hand 
the conductance properties of interacting quantum dots 
with a superconducting lead have been discussec&i 9 - - — . 
On the other hand the full counting statistics (FCS) of 
multi-terminal systems with a chaotic cavity and a su- 
perconducting connector were investigated^ 2 - - — . Among 
other things the cross-correlation noise has been found 
to be an efficient tool for discrimination of the current 
caused by CAR from other current contributions^—. 



taken into account. In this paper we would like to close 
this gap by considering a quantum dot modelled by a 
noninteracting resonant level. This setup enables one 
to access the full energy dependence of transmission 
characteristics generated by the superconducting density 
of states (DOS) as well as by the energy-dependent 
transmission of the dot. 

The quantity of our primary interest is the FCS in 
terms of the cumulant generating function (CGF) \- 
It represents a very convenient tool for the calculation 
of a variety of transport properties. It is directly re- 
lated to the probability distribution P{Q) to transfer Q 
elementary charges during a fixed very long measure- 
ment time r. By a simple derivation with respect to 
some parameters (counting fields) x gives all cumulants 
(irreducible momenta) of -P(Q) 36 ' 37 - The relevance of 
higher order cumulants has been recently demonstrated 
experimentall y 38 : 39 . Moreover, the analytic structure of 
the CGF provides deep insights into the different trans- 
port processes^— . As we shall demonstrate later, in 
the case of the multi-terminal quantum dot considered 
here, the FCS will reveal the different facets of entan- 
glement that we expect to be observable in upcoming 
experiments. 

The Hamiltonian for the system under consideration is 
given by 



H — Hq + Hd + Ht + H s 



(1) 



Even in its basic realization the noninteracting quan- 
tum dot usually shows a highly nonlinear transmission 
leading to nonlinear current- voltage relations. Combined 
with a superconducting electrode, the single-particle 
density of states of which is itself highly nonlinear due 
to the finite gap, the resulting structure is expected to 
possess transport properties with a nontrivial voltage 
dependence. Thus far these aspects have not been fully 



The term Hq describes normal or ferromagnetic elec- 
trodes in the language of the respective electron field 
operators ^ aa {x), where a = 1, . . . , N numbers the elec- 
trodes. The local density of states po a in the leads is 
assumed to be very weakly energy dependent in the rel- 
evant range of energies. The normal electrodes are mod- 
elled by noninteracting fermionic continua held at the 
chemical potentials jj, a . Ferromagnetic electrodes are de- 
scribed by the Stoner model with an exchange energy h cx 
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Consequently they can be described as fermionic con- 
tinua with a spin-dependent DOS pa aa = Po«(l + C-Pa), 
where P a is the polarization 45 . The superconducting 
electrode is described by the BCS Hamiltonian with the 
gap A of the superconducting terminal 4 ^ 



H, 
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k.a 



The superconductor is kept in equilibrium as in previous 
treatments of similar problems 41 *. The applied bias volt- 
ages are given by V a = fi s — (j, a = —fj, a (we use units 
where e = h = ks = 1)- The electron exchange between 
the dot and the electrodes is given by 4 ^ 

H T = E^ K^aaix = 0) + h.C.] 

a,(T 

+ ^ 7s [4* SCT (x = 0)+/ l .c], 

where "/ a or 7 S are the tunneling amplitudes between the 
dot and the normal or superconducting electrodes. The 
tunneling is assumed to be local and to occur at x — 
in the coordinate system of the respective electrode. d a 
is the annihilation operator of an electron with spin a on 



the dot. The quantum dot in the presence of a magnetic 
field B is modelled by H d = £ CT (A d + crh/2)d+d a =: 
Scr dvdtda, where is the bare dot energy and in SI- 
units h — fJ-sgB with Bohr's magneton fiB arid the gy- 
romagnetic ratio g^. 

In order to determine the FCS we calculate the CGF 
lnx(A) = ln(e iA< 2), which depends on the counting fields 
A = (Ai, . . . , Ajv, A s ) in the respective electrodes. It pro- 
vides all the higher statistical moments (or irreducible 
cumulants) «Q n )) of P{Q). The CGF is calculated us- 
ing the Keldysh Green's function approac h 36 i 49 adapted 
to quantum impurity problems in^2r— . According to the 
generalized Keldysh approach one obtains for the CGF 40 



X(A) = (T c exp 



dm, 



A(t) 



C 



where the dependence on the counting fields is contained 
in 



Hr [t) = J2^[ eiXa(t)/2d ^ a ^ X = ° ) + h - C - 



E%[< 



iX *W 2 dt*s<r(x = 0) + h.C. 



The different counting fields are nonzero only during 
the measurement time r and have different signs on 
the forward and backward path of the Keldysh contour. 
The scattering matrix has two energy regimes^ with 
energy-dependent transmission coefficients. The respec- 
tive Fermi distributions of the individual terminals are 
abbreviated by n a and n a+ := 1 — n a (— uj) for hole-like 
contributions. This allows to express the CGF in the 
form given in equation (J2]). 



lnx(A) = — I duj 

7T J 



U-A 



E 111 \ 1 + E T ij<T (oj) ni (l - njXe^-^ - 1) 

i,j=l,...,N,s, 



1„ / A- |.x 
— i 
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E ln ' 



!+ E TAijreTliil -Uj)^-^ -1) 

i,j=l,...,N, i=£j 



1 + Yl T AijtTh n j+ (l - n t+ )(e l ^- x ^ - 1) 

i,j=l,...,N, ij&j 



T Aia \ ni (l-n i+ )(e 2i ^-^-l) 



=1,— ,iV 



+ X) TcM ^ h+( X - n t )(e^-^-^ - 1) + m(l - n j+ )(e-« 2X >- Xi -^ - 1) 

i,j=l,...,JV, ijtj 



, (2) 



where we define 



r y<r (w) = 4r J r J (i + ap){i + aPj)/[(" - 5 °f + ( E Tk ^ + aPk ) + r si) 2 ]' with ( p * = °) ( 3 ) 

k=l,...,N 

T Aijae = ATdl + aP i )T j (l + oP j ){{lo-5^) 2 + [ £ T k {\ - aP k )] 2 } / det Aa (u), (4) 

fe=l,...,JV 
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T Aijah = 4^(1-^)^(1 - oPj){(u - 8 a ) 2 + [ ^ r fe (l + CT P fe )] 2 }/deW(w), (5) 

k=l,...,jV 

T Aia = 4T 2 (1 - aPi)(l + CTPi)rydeWM, T C ^ CT = 41^(1 + aP^T^l - aPj)r 2 s2 /det ArT (uj), (6) 
deWM = (w - <5 CT ) 2 (w - 5- a f + {[ ^ T k (l + aP k )} 2 +T 2 s2 } 

fc=l,...,JV 

x{[ £ iMi- CT p fe )] 2 + r 2 2 } + [ £ r^i + ^V-^) 2 

fc=l,...,jV fc=l,...,jV 

+ [ £ (l-aP,)] 2 (c- ( 5 ff ) 2 + 2r 2 2 ( W - ( 5 CT )( W -5_ ff ), (7) 

I 



k=l,....N 



The abbreviation = 7rpoi|7i| 2 /2 is the (energy- 
independent) dot-lead contact transparency with dimen- 
sion energy for the normal leads. For the superconduct- 
ing leads it is affected by the energy- dependent super- 
conducting DOS so that 



given by 



r sX = irp as \"f s \ 2 \u>\/ (2\/oj 2 - A 2 ) 
r s2 = 7rp 0s | 7s | 2 A/(2VA 2 -^). 



(8) 
(9) 



For \u)\ > A the counting factors e^ Xi ~ x ^ describe single 
electron transport between the different terminals. For 
\lu\ < A the superconducting DOS only allows excita- 
tions from the Cooper pair condensate which leads to 
different transport characteristics. They are described 
by counting factors e 2l ( A s- A ;) anc } e 24A s -iAi-jAj re f err i n g 

to the transfer of two particles from the superconductor 
to a single or two separate terminals. The transmission 
coefficients T A ia and Tc A ij a therefore refer to Andreev- 
reflection and CAR respectively. 

For A — > the CGF reduces to the well-known expression 
for the CGF of a noninteracting dot in a multi-terminal 
geometry^. The same holds for js —> 0. 
The FCS for a three-terminal structure at T = us- 
ing a chaotic cavity with energy-independent transmis- 
sion instead of a resonant level have been calculated 
befor o 22 ' 26 . In this case the CGF adopts a characteristic 
double square root form. The first square root instead 
of a logarithm is due to the diffusive transport through 
a chaotic cavity^. The second square root may be ex- 
plained looking at physical observables that are calcu- 
lated via derivatives of the CGF, where it leads to addi- 
tional factors 1/2. Therefore it corresponds to the fac- 
tors 1/2 in the transparencies Ti, T s \ and r s2 that are 
due to the separate treatment of electrons and holes as 
a consequence of the proximity effect. The transmission 
coefficients arc different in our case because of the energy- 
dependent DOS of the quantum dot. As in our case the 
ones for single electron transmission, direct Andreev re- 
flection and crossed Andreev reflection are proportional 
to rir 2 , (r 2 +T 2 )r 2 and 2rir 2 r 2 , respectively. 
Considering the case of only a single normal electrode 
the conductivity at low voltages and no magnetic field is 
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In the last step we introduced fj = 2Fi, f s = 2T S and 
restored Si-units. The result coincides with the one pre- 
viously obtained in££. 

The three-terminal case with two normal metal drains is 
of particular interest. Positive cross correlations between 
the normal electrodes can be used to probe the existence 
of entanglement^. Along with the noise the cross cor- 
relation can be calculated as a second derivative of the 
CGF 



P 1 



1 d 2 ln X (A,r) 
t d\\d\2 



A=0 



Depending on the coupling and voltages choice differ- 
ent cross correlations may be observed. Three differ- 
ent types of transport between the normal drains are 
present: the direct single electron current proportional 
to rir 2 (Vi - V 2 ), the direct Andreev reflections (DA) 
proportional to 2(r 2 Vi + r^V^r 2 or CAR proportional 
to 2Fir 2 r 2 (Vi + V2). If the superconducting terminal 
is weakly coupled to the quantum dot, cross correlations 
will either be dominated by single electron transmission 
or, for Vi ~ V%, by direct Andreev reflections to a nor- 
mal drain leading to negative correlations. There are two 
cases in which positive cumulants are observed. If the su- 
perconductor is coupled better to the quantum dot than 
the normal terminals, positive correlations may be ob- 
served for voltages close to the superconducting gap and 
V\ w — V 2 , see Figure [TJ In this case CAR is strongly 
suppressed and one expects single electron transmission 
to be dominant. However, the energy-dependent DOS 
of the superconductor leads to large transmission coef- 
ficients for double Andreev reflections from one normal 
drain to the superconductor and further to the second 
normal electrode. 

This Andreev-reflection Enhanced Transmission (AET) 
is also observed if an additional broadening of the BCS 
DOS is taken into account as mi 2 -, where ui — > oj+irj s with 
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FIG. 1: Cross Correlations P/ 2 calculated from the CGF in 
equation ((2| with the parameters T s = A/2 = Ti, r 2 — 
0.2A, A d = 0, T — 0.1A as a function of V 2 for Vi = -V2 
(solid line), Vi = (dotted line) and Vi = V2 (dashed line). 
The three cuts show that one observes negative cross corre- 
lation for Vi = V2 and Vi = whereas a positive cross corre- 
lation maximum for Vi = —Vz close to the superconducting 
gap caused by AET is obtained. 



a typical experimental value r\ a of about 1CP 2 A. AET is 
a robust phenomenon. In modern experiments using ei- 
ther InAs nanowires- or carbon nanotubesZ T s rs A is 
generically obtained. The hybridisation with the normal 
leads can be tuned via top gates to e.g. T n T s /2. AET 
should then be observable via cross correlations or from 
the direct currents. 

The second case for positive correlations is observed 
for strong coupling of the superconductor to the quan- 
tum dot and asymmetrically coupled normal terminals. 
Since the contribution by CAR is proportional to V\ + V% 
a voltage bias may be applied via the the weakly coupled 
normal electrode j — 1,2. The DA contribution is pro- 
portional to Tj and therefore will be strongly suppressed 
for this weakly coupled drain. For voltages well below 
the gap CAR dominates over single electron transmis- 
sion because a CAR is possible in two ways. In Figured 
we show the cross correlation in this situation. 
The directions of the currents are different in the case 
of CAR and AET. CAR describes electron transfer 
into/from the superconductor whereas in the case of AET 
the superconductor only assists electron transfer between 
the normal drains. The FCS formalism allows us to follow 
the different transport processes independently. In the 
case of CAR two electrons from the same Cooper pair 
are transferred to spatially separate electrodes. In the 
case of AET an electron impinging on the superconduc- 
tor is retroreflected as a hole. The same Cooper pair that 
was generated by this Andreev reflection is then trans- 
ferred further by a hole from the second normal drain 
that is retroreflected as an electron. In the case of CAR 
we therefore observe entangled electrons in the normal 
electrodes both in energy and spin space. Considering 
AET the electron and the hole originate from one Cooper 
pair and due to the conservation of spin and energy must 
also be entangled in energy and spin space. Only AET 
and CAR generate positive correlations which therefore 
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FIG. 2: Cross correlation P/ 2 calculated from the CGF in 
equation @ with the parameters F s = A, Y\ = 0.4A, A^ = 
0, T = 0.1 A, Vi = as a function of V 2 . The dotted line 
shows the result for r 2 = 0.05 A, the solid line is for r 2 = 0.1 A 
and the dashed line is for F 2 = 0.15A. Consequently the sec- 
ond electrode is weakly coupled. Therefore the CAR processes 
dominate for voltages below the gap over single electron trans- 
mission and DA leading to a positive cross correlation. This 
effect is enhanced for weaker coupling of the second electrode. 
The effect is also present at finite Vi but it is weakened. 



indicate the presence of entanglement. 
The asymmetric coupling necessary to observe CAR can 
easily be realized experimentally by fabricating a simple 
NQS structure on a substrate with a single normal drain 
and using an STM tip as the second weakly coupled ter- 
minal. 

There are two further possibilities to observe CAR as 
the dominant transport channel. On the one hand us- 
ing additional normal drains enhances the possibilities 
for CAR as can be seen directly from the CGF in equa- 
tion ©. The combinatorical factor naturally leads to a 
dominance of CAR. On the other hand DA and single 
electron transmission can be fully suppressed by ferro- 
magnetic drains with antiparallel polarisation^ 5 -. In this 
case we always observe positive correlations for voltages 
below the gap and not too high temperatures if the po- 
larisations are chosen strong enough. For P\ = — P2 = 1 
we also considered the case of three normal terminals 
and we obtained P/ 2 — in accordance with previous 
calculations^ 6 -. 

A ferromagnet preferably accepts electrons of a specific 
spin. It may therefore also be viewed as a detector for 
this specific spin directio n 25 ' 56 . This fact enables us 
to demonstrate a possible violation of the Bell-Clauser- 
Horne-Shimony-Holt inequality along the lines o f 25 ' 56 . 
First we have to include the spatial direction of mag- 
netization into our previous treatment, where we only 
considered one specific quantization axis for the polar- 
izations of the ferromagnet. 

We consider a device as in Figure [3] with four F n (n — 
1, • • • ,4) terminals. The drains Fi,Fq, and ^3,^4 are 
pairwise polarized in opposite directions to serve as spin 
detectors. First we only consider the case with a sin- 
gle quantum dot in the middle. For a simplified discus- 
sion we want to consider the special case P\ — —Pi — 
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P 3 = -P 4 and Tx = T 2 , T 3 = T 4 as in 56 . To avoid 
confusion with the different quantization axes we limit 
our discussion to h = 0. The energy- dependent trans- 
mission coefficients for CAR in the case of four ferro- 
magnetic drains and a — 1,2, /3 = 3,4 are given by 
TcA a /3* = 4r x r 3 (l - PiP 3 )rydeW(w). The polar- 
ization in the direction-dependent case may be rewrit- 
ten as Pi — >• gi = Pim.i, i = 1,3 with a unit vec- 
tor m.i describing the magnetization direction^. In this 
way we obtain the transmission coefficients TcAufio = 
4Tir 3 (l - g Q g^)r^ 2 /det J 4 cr (oj). The differential conduc- 
tances at T = for the respective crossed Andreev re- 
flections to drain a and /3 can be directly read off via 
G C A a /3a = & [T CA *pa{V)+Tc Aa p*(-V)]. As CAR is 
a nonlocal process these conductances describe nonlocal 
correlations. 




FIG. 3: Schematic setup for demonstrating the violation of 
Bell's inequality. A nanowire is connected to four ferromag- 
netic drains fi,...,4 and a superconductor S. The top gates 
TGi,2 may be used to divide the nanowire into three quantum 
dots QDi,2,3 with tunable coupling. 

They also give the probabilities p a fi for simultaneous de- 
tection of an electron at a and /3 by 

Va.fi = GcAa(3a/ ^ GcA-fSa ■ (10) 
{7,5} 

The drains 1,2 and 3,4 represent spin detectors (Alice, 
Bob) with respect to two possible choices for the direction 
gi,3i g'i3- The spins can be detected by the respective 
currents. We shall discard events with both electrons go- 
ing to the same detector by normalizing the probabilities 
to go to different detectors as in££. 

We should stress that we need to record every single event 
separately. To go over to a non-time-resolved measure- 
ment scheme one needs to observe not just the spin cur- 
rent but the spin-current fluctuations^. 
As i n 25 ' 56 we find for the Bell parameter e = |gig3 + 
g'ig3 + g'ig3 — gig3 1 ■ The maximum of e is 2^/2. Bell's in- 
equality e < 2 is violated if the polarization ||gi|| = Pi = 
P 3 > 2~3 as in other treatments^. Such a degree of po- 
larisation is hard to obtain with ferromagnets but may 
easily be reached using double-quantum dot structures— 
This also represents a direct measurement of the concur- 
rence (C) since the maximal value of the Bell parameter 
tmax = 2vl + C 2 , where C is a measure of entanglement 



with possible values between and 1. 

The same expression for the Bell parameter has been 
obtained previously by Morten et al.— for a chaotic cav- 
ity. Our result shows that the concurrence does not de- 
pend on how the beam splitter geometry is realized. The 
cross correlation however is affected by the difference in 
the transmission coefficients and therefore only indicates 
the presence of entanglement. The described experiment 
appears to be a perfect analogue to the optical experi- 
ments using parametric downconversior^. However, care 
has to be taken with respect to the distinguishability of 
particles. In optical experiments one detects spatially 
separated particles where entanglement is unambiguously 
defined. In the case of a single quantum dot we detect 
the spin direction of electrons coming from the dot which 
are indistinguishable. In this situation the state has the 
form of a single Slater determinant, which is considered 
a nonentangled stated. But we can easily go over to 
the distinguishable case by using the full power of the 
experimental setup in Figured! By tuning the top gate 
voltages one can induce tunnel barriers to obtain a triple 
quantum dot in series^. If we go over to a strongly cou- 
pled triple dot we obtain well separated densities of states 
and therefore detect distinguishable particles from QD\ 
or QD 3 . In this case we can neglect resonant tunneling 
and evaluate the transmission coefficients by multiplying 
the transmissions through the different quantum dots. 
In this situation we also obtain the same Bell parame- 
ter, however now for an unambiguously entangled state. 
The correlations remain the same in the distinguishable 
or indistinguishable case. 

To conclude we have analytically calculated the full 
counting statistics (FCS) of a noninteracting quantum 
dot (resonant level) contacted by multiple metallic elec- 
trodes in the normal state as well as by a superconducting 
one. The superconducting ground state leads to posi- 
tive cross correlations between the currents in the nor- 
mal drains by two transport phenomena that we identi- 
fied as crossed Andreev reflections (CAR) and Andreev- 
reflection enhanced transmission (AET). Our results al- 
low for an analysis of a possible violation of the Bell- 
Clauser-Horne-Shimony-Holt inequality that is related to 
a measurement of the concurrence C. Interestingly, we 
find it to be independent of the constellation of coupling 
parameters. Moreover, it coincides with the results found 
for the quantum dot realization via a chaotic cavity. Fi- 
nally we discussed the question of entanglement for in- 
distinguishable and distinguishable particles. We have 
shown that the correlations should be unaffected by this 
distinction and how both cases can be addressed experi- 
mentally. 

The authors would like to thank S. Maier and P. Plotz 
for many interesting discussions. The financial support 
was provided by the DFG under grant No. 2235/3, 
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